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Traffic delay, caused by temporary all-lines-busy conditions, is analyzed 
for three mathematical switching models. They are classified as "address 
camp-on," '^retrial," and "message storage" models. The models are de- 
signed to permit a study of basic traffic theoretical problems encountered in 
the rapidly growing field of data communicatio7is, but they are not identical 
with any of the existing data switching systems. Each model assumes that 
a message is switched only ihroiigh one switching center which must es- 
tablish connections via line groups to one or more addressed receiving sta- 
tions, i.e., each model contains only a single switching center. Numerical 
residts for the average delay on all messages are obtained on the IBM 7090 
computer. 

I. INTRODUCTION 

Switching centers can be used to link together communication lines 
for the transmission of data between a variety of business macliines and 
computers. Due to randomness in the required interconnections a switch- 
ing center may occasionally not find an idle line to a particular receiving 
station, so that a delay can occur. More than one method can be followed 
when a switching center finds all lines to a receiving station busy. Some 
switching models appear to obtain lines to the addressed receiving sta- 
tions in a shorter time than other switching models. This means that 
with one switcliing model a gi\'en delay requirement can be met with 
fewer lines than with another switching model. Tliis is not to say that 
the model which would render a given grade of service witli the least 
number of lines also is the most desirable from an economic point of 
view, because delay is only one factor which enters into the choice be- 
tween data switching systems. Components of a switching system, such 
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as lines, memory, etc. do not bear the same price tag, and minimizing the 
number of components of one kind does not ensure economical efficiency. 

Interest in the particular traffic engineering problems of data switch- 
ing has been present for at least 15 years. Yet, traffic engmeering work 
was mainly concentrated on classical telephone trunking problems, and 
a variety of such fundamental problems have been worked out. Some 
of the data switching systems which are being studied or are now in use 
cannot be analyzed by standard mathematical approaches of traffic 
theory because the operating conditions differ from those of the mathe- 
matical models used in the analysis of classical telephone problems. The 
understanding of the fundamental traffic theoretical problems encoun- 
tered in data switching is a prerequisite for an exact mathematical analy- 
sis of message delay in data systems. The fundamental problems need 
to be studied on simplified models which lend themselves best to mathe- 
matical treatment and, therefore, will not be identical with any of the 
present data systems in use. We have constructed for study three hypo- 
thetical models which we cafi "address camp-on," "retrial," and "mes- 
sage storage" models, and have analyzed message delay for each of them. 

Message delay is defined as the delay between initial request by the 
switching center for a line, and the moment the message is released 
from the switching center for transmission. The switching center handles 
messages in a manner described by one of the three switching models. 
The delay is caused by temporary all-lines-busy conditions in the line 
groups which connect the switching center with the addressed receiving 
stations, This type of traffic delay must not be confused with the total 
delay from the time a message is ready at the data source and the time 
the message is actually received at a destination. No account is taken 
of messages which arc switched through more than one switching center 
in tandem. 

This study, then, shall not be looked upon as an attempt to make a 
choice between switching systems, since such a choice cannot be based 
solely on the delay performance of mathematical models. A true com- 
parison between switching systems must include other factors, as for 
instance the cost of memory and logic, loading of transmitters, and load- 
ing of incoming fines, all of which are neglected here. 

II. DESCRIPTION OF MATHEMATICAL SWITCHING MODELS 

The following describes each of the three switching models. The de- 
scription is preceded by an outline of features which are conunon to 
each model. The mathematical derivations given in the appendices and 
the delay curves are based on these models. 
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2.1 Common Features 

Think of a data source feeding messages into a smtching center that 
has a large number of line groups radiating from it (see Fig. 1). 

Each line group connects one, and only one, recei\'ing station with the 
switching center. A receiving station is capable of receiving from all 
lines in its group simultaneously. There might be one or more lines per 
group, but each group contains the same number of lines c. Full access 
is given to each line in a group. A message is said to have .1 addresses 
when a copy of the message must be transmitted over A different line 
groups to A different receiving stations. The number of addresses per 
message remains constant for all messages. The switching center is 
responsible for transmitting a copy of the message to each of the ad- 
dressed receiving stations. The addresses of a message are cho.scn at 
random from a large number of possible receiving stations. This permits 
us to assume that all line groups are independent of each other. Mes- 
sages are originated and addressed in such a way that Ui , the information 
load offered to a group, is the same for every group. The information 
load is defined as the number of first, i.e., unrepeated, message attempts 
which are expected to be generated during an interval equal to one aver- 
age message length. First attempts are made Poisson distributed hi 
time, meaning that the probability that exactly k first attempts are 
generated during an interval of length I is given by 
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Fig. 1 — Switching model. 
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Pr (k,t) = 



k\ 



in which t is in units of the average message length. The length of all 
messages is exponentially distributed with mean 1. In the mathematical 
derivations the average message length is taken as the miit of time. An 
exponential message length distribution is chosen, since it is believed 
that it will serve as a good approximation in a larger number of practical 
cases than a constant length. The instant a new or repeated message is 
originated, the switching mechanism begins to hunt for an idle line to 
each of the addressed receiving stations. No delay is imposed by the 
switching mechanism itself. Each new message is eventually delivered to 
the respective receiving stations, i.e., no messages are lost. The system is 
in statistical equilibrium, which is to say that the system is in the steady 
state such that the average number of messages in the system during 
any long interval of time remains constant. 

When all lines in a group to one or more of the receiving stations are 
busy wc nnist find a way of delaying delivery. The camp-on and storage 
models assume that blocked requests form queues at the switching cen- 
ter. The retrial model assumes that blocked re([ucsts are withdrawn from 
the switching center and reofferod at a later time. 

The reader has doubtlessly observed the very simplified and idealized 
set of common features on which the switching models are based to per- 
mit mathematical analysis. The same applies to the features which are 
unique to each of the three models. 

2.2 Address Camp-On Model 

When a line group to an addressed receiving station is blocked, the 
request for service in this group will camp-on and wait in the order of 
arrival until a line is assigned by the mechanism which scans continu- 
ously for idle hues. The assignment of available lines to waiting requests 
is done on a "first come, first served" basis. When a line is assigned it is 
immediately made busy. The message, however, is not released from the 
data source until lines to all addressed receiving stations are secured. 
When the last line is obtained, the message content is released and trans- 
mitted simultaneously to each of the addressed receiving stations, after 
which the fines are released. The holding time of a line in the camp-on 
model is made up of the sum of two random variables: namely, the ex- 
ponentially distributed message length and the time spent waiting vmtil 
lines are secured to all addressed receiving stations. 
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2.3 Retrial Model 

In the retrial model, a message is released from the data source only 
when linos are found at the switching center to all addressed receiving 
stations. If even one line group is blocked, the message is not transmitted 
to any of the addressed receiving stations and is temporarily cleared from 
the switching mechanism without making any lines busy. A blocked 
message is reoffered any number of iin.es from the data source after a 
constant time interval t, luitil an idle line is found simultaneously to 
each addressed receiving station. At a retrial of a blocked message, an 
attempt to seize an idle line is made in the same groups as at the previous 
attempt. The message delay is determined by the number of attempts 
made and by the length of the constant retrial interval r. The holding 
time of a line in the retrial model is equal to the exponentially distributed 
mes.sage length. 

Another way of making retrials is to let the delay in the delivery of 
the message content to any one addressed station be independent of the 
d(^lay to the other addresses of a multiaddress message. In this case, a 
message having A addresses would be considered to consist of .4 one- 
address messages and the delay would be that given for the one-address 
case of tlie retrial model. 

2.4 Message Storage Model 

Message storage is analysed on a model in which requests for lines in 
a busy group form (|ueues in the order of arrival. In the multiaddress 
case, some addresses of tlie message may fmd their line groups busy 
while other addresses may obtain lines to the addressed stations with 
no delay. The model assumes that the message is released with no delay 
to stations which are not blocked, and that the delivery to a blocked 
station is delayed only until the instant a line is found by the switching 
mechanism which scans continuously for idle lines. As in the retrial 
model, the line holding time is eciual to the exponentially distributed 
message length. 

III. METHOD OF ANALYSIS 

Tlie delay performance is analysed as messages are switched through 
one switching center which employs one of the three described switching 
models. It should be pointed out that the results obtained here apply 
only to the niathematical models used. All approximations mentioned 
in the analysis are approximations of the model to which they refer. 
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The mathematical analysis of the delay performance of the address 
camp-on model is given in Appendix A. The problem is to find the total 
occupancy of the outgoing line groups. In the multiaddress case, out- 
going lines can be held busy in excess of the message length. This excess 
holdmg time mcreases the load carried from the useful information load 
ai to a total load Or . The excess holding time is the average length of 
time between line seizure and the time lines are found for all addresses 
of the message. Erlang delay probabiUty is used by the introduction of 
an approximation which assumes that the total holding time of an out- 
going line is exponentially distributed. No explicit expression is derived 
for the total load ar . Solutions for Or are found by solving (7) and (10) 
of Appendix A in an iterative computer program. The average delay 
follows from (11). 

In Appendix B the mathematical analysis is given for the retrial 
model. The retrial method has been imder consideration for application 
in both mihtary and commercial data systems, and this method is also 
used in voice telephone communications. The mathematical analysis of 
delay in systems in which blocked attempts are reoffered is one of the 
fundamental traffic problems for which an exact solution is not available. 
The prospect of using the retrial method in data systems emphasizes the 
need to treat such systems analytically. The analysis given here is 
not exact because a number of approximations had to be mtroduecd 
to obtain numerical results. Since the retrial method is a basically 
unsolved problem it must first be studied in its sunplest form, which 
exists for the case of one address per message. Considerable effort, there- 
fore, is spent in Appendix B on the discussion of the one-address case. 
Our approach to the retrial problem is to find approximations for the 
unconditional state probability of finding i Imes in a group of c lines 
busy, ^ z ^ c. Then, approximations are foimd for the conditional 
probabilities of finding i Imes busy at ^o + r, when the state of the group 
is known at U , k - r, U - 2r, etc. The delay for the one-address case 
follo^vs from (23) of Appendix B. For the three-address case the delay 
is computed from (31) in a Markov process which is m itself an approxi- 
mation of the retrial problem since it accoimts only for a first-order 
dependency. 

The basic problem m the retrial model is to find approximations of the 
conditional probabilities mentioned above. These are obtained by in- 
tegratmg a set of differential equations (16), using a line request rate 
o}(t) which by itself is conditioned on previous states of the line group 
and, therefore, is dependent on tune. The line request rate a)(i) appears 
as a coefficient m (16). Smce wit) can be expressed only as a function of 
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solutions to (16), we cannot find o}{t) explicitly, but must compute it in a 
long process of progressive iterations. It will become apparent from Ap- 
pendix B that not all appro ximationf; made can l)e dearly justified, but 
the results olstained are sufficiently accurate for comparison with other 
switching models. Some of the approximations appeal- ci'itical for short 
retrial intervals t, particularly when c is small. The amount of effort 
and computer time spent on solving the retrial problem analytically is 
not necessarily less than the amount of effort and time spent by simula- 
tion. The problem is l)y no means solved, but it is hoped that by this 
analytical approach the way is paved toward a more complete analysis 
of retrial systems. 

For message storage, the average delay can be determined by the well- 
established methods of traffic theory developed by A. K. Erlang. These 
are outlined in Appendix C. The a\'erage delay for the storage model is 
computed from (32) of Appendix C, and no approximations need to be 
made. 

IV. RESULTS AND CONCLUSIONS 

For a fixed amount of information load, each switching model pro- 
duces different delays. This means that some switching models must be 
operated at lower occupancy than others to ensure tlmt delays encoun- 
tered will not exceed the desired maximum. The delays sliown below 
for each switching model do not necessarily keep their relationships in 
respect to each other wlien messages are switched tlirough several switch- 
ing centers in tandem. 

The results of computations for one address per message arc shown m 
Figs. 2 and 8 for one and ten lines per group, I'cspectively. I'igs. 4 and 5 
show similar results for three addresses per mesi^age and one and ten 
lines per group, respectively. The "information occupancy" in these 
figui-es is numerically equal to the information load olfered to the line 
group divided by the number of lines per group, i.e., ai/c. The term 
"occupancy" refei's to the percentage of time a line is occupied on the 
average. The fraction of time a line is actually utilized for the transmis- 
sion of information, then, is er^ual to ai/c, so that we may also call "in- 
formation occupancy" the "line utilization." 

First let us discuss the address camp-on model. This method offere 
the advantage that eiror correction can be performed on multihnk con- 
nections on an end-to-end basis because the message content remains in 
storage at the data source until a connection is set up to all addressed 
receiving stations. The camp-on model also is of interest because storage 
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0.3 0.4 0.5 0.6 0.7 

INFORMATION OCCUPANCY 



Fig. 2 — Average delay on all messages vs information occupancy; 1 line per 
group, 1 address per message. 

at an intermediate switching center need be provided only for the ad- 
dress portion of a blocked message and this might have some economic 
advantages over other switching models. For one address per message 
there is, in theory, no difference in the traffic dohxy performance of the 
address camp-on and me.ssagc storage models when the message is 
switched only once. It should be remembered that in the camp-on mode! 
an intermediate switching center keeps the incoming lines busy in excess 
of the message length for the duration of a delay which, for a given 
information load ai, increases the actual load carried. Because we con- 
sider only single-switched messages, no account is taken here of this 
type of line loading. 

As was mentioned before, the total holding time of an outgoing line 
in the address camp-on model is made np by the excess holding time, 
whicli is the time spent waiting for other addresses to find lines, and by 
the actual message length. In Fig. 6 we show the total occupancy ar/c 
versus the actual information occupancy or line utilization Or'c for three 
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addresses per message. We sec that the total occupancy approaches 100 
per cent at a surpiisingly low information load. This is due to the fact 
that the excess hfilding time increases the load on the outgoing line 
groups, which in turn increases delays and thus brings about longer 
excess holding times. This makes the camp-on model unusable beyond 
certain intolerably low le\'els of line utilization. For instance, in the 
three-address case, line utilization must be limited to about 14 per cent 
or 00 per cent for line group sizes of c = 1 or c ^ 10, respectively. It 
can be seen in Fig. fi that beyond this point the total occupancy blows 
up and with it the delay imposed on a message. A similar result was 
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Fig. 4 — Average delay on all mesaagos vb inrormalion occupancy; 1 line per 
group, 3 addressea per message. 

obtained by Weber^ in a different approach. A simulation made by 
Weber' for c = 1 shows close agreement with the results derived here. 
For example, the simulation shows a maximum utilization of about 14 
per cent for c - 1 and three addresses per message, which is the same 
as derived here analytically. This indicates that considerable confidence 
may be placed in the approach presented in Appendbc A. 

From the delay performance of the address camp-on model it is con- 
ckided that any switching method in which delays become a substantial 
part of the Ime holding time will require a relatively large number of 
lines to provide adequate service. By the same token we may conclude 
that even more lines will be required when the message is switched more 
than once, i.e., through more than one switching center. 

Next we turn our attention to the retrial model. The performance of 
the retrial model as a function of the retrial interval, r, is of interest. 
We observe in Figs. 2-5 that when one doubles t the delay is less than 
doubled. On the one hand, we expect longer retrial intervals to cause 
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longer delays, On the other hand, it can be shown that the probability 
that the message succeeds on a retrial increases with increasing length 
of T. vShorter retrial intervals result in smaller chance for success than 
longer retrial intervals, but the fact that in any given time there arc 
more attempts made with short retrials than with long reti-jals makes 
the average delay a monotone increasing function of r. For large values 
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pjg 6 _ Xotiil occupancy vs information occupancy for address camp-on 
model; 3 addresses per message. 

of T, say at least twice the message length, the rate of increase of the 
delay will be almost proportional to the rate of increase in t. This is 
so because the rate of change of the probability of being blocked again 
becomes smaller with increasing values of r. With t approaching 0, the 
a^'erage delay with the retrial model should approach the delay for the 
message storage model in the case of one address per message. 

The delay in the retrial model increases with the number of addresses 
per me.ssage. A comparison between Figs. 2 and 4 shows that this in- 
crease is quite significant when the number of lines per group is small. 
This increase can, of course, be avoided when a multiaddress message is 
broken up into se\'eral one-address messages, as suggested earlier in Sec- 
tion 2.3. Stations not blocked would then receive the message content 
independent of the availability of lines to the other stations. If this mode 
of operation is used for retrials, the retrial curves given in Figs. 4 and 5 
are to be replaced by those for the one-address case shown in Figs. 2 and 
3, respectively. Such a change of the retrial model would bring the av- 
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orage delay lor \'ery short retrial intervals close to that oi' the message 
storage model, independent of the number of addresses. 

The retrial model has the advantage that no sf.orage hats (.0 be pro- 
vided in a switc^hing center except at the data source. However, provision 
must be made to instruct the data source to reoffer the message when an 
all-lines-busy condition is eneoimtered at a switcliing center. The retrial 
method is particularly well suited for error correction by retransmission 
from the data source on request from a terminating station which de- 
tected an error. 

Since the analysis of the retrial model is based on approximations, we 
do not expect absolute accuracy of the curves derived. The retrial model 
has been simulated by others for some special cases and it is interesting 
to compare the results. This is done in Table I, in which time is expressed 
in units of the average message length. We observe some disagreement 
between analysis and simulation for large retrial intervals. 

Finally, let iis look at the message storage model. Tliis method pro- 
vides an efficient use of lines, even if the line groups are small, and delay 
is independent of the number of addresses per message. It requires, how- 
ever, that considerable storage be provided because each switching 
center must be arranged to permit full message storage to allow for de- 
lays exceeding a message length. Provision must also be made for trans- 
mission of a copy of the message to each addressed station independent 
of the transmission to any other station. 

According to the model of the message storage method, a message is 
delayed only imtil the very instant a line is found. From a practical 
point of view this means that the line-hunting mechanism should be 
activated as soon as the address is decoded. On the other hand, message 
storage may be operated so that the message is stored completely at the 
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switching center before the line hunting starts. In this case, the amount 
of time needed for full message storage must be added to the delay. The 
latter mode of operation adversely affects the delay performance of the 
storage model, particularly when the message is long in comparison with 
the delay that can be tolerated. 

The delay performance for the storage model will become consider- 
ably worse than shown when the line back to the originating station 
cannot be released as soon as the message has been transmitted over it. 
For multiswitched messages the release of lines between switching cen- 
ters would ensure that the line holding time is not increased by the de- 
lay. 

The curves given for the storage model can be considered accurate 
because the vahdity of Erlang delay formulas has long been observed. 
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APPENDIX A. 

Mathematical Analysis of Address Camp-On Model 

Let the excess holding time Ea be the average time between seizure 
of a line and the time lines are found for A addresses of a given message. 
Further, let Zk,± be defined as the average time between initial request 
for lines and the time K out of A addresses have seized lines. At the time 
A addresses have seized hnes, the message is ready to be transmitted. 
The average delay on all messages, d^ , is defined as the average time 
between initial request and the time lines have been seized by all A 
addresses, as illustrated in Fig. 7. 

The average is a hnear operator, and one obtains for the expected 
excess holding time 

A K=l 

And since 

Wk,a = dx — Zk.a 
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it follows that 



Ea = d,-\'tz^.,. 



(1) 



1 "^ 
The term -r "Y^Zk.a is recognized as the average time between re- 

A K^l 

quest for and seizure of a line for any given single address. 

We define Q{t) as the probability that the delay is less than or equal 
to / between the time of a request for a line by a given single address 
and the time of line seiztu-e, and obtain 



A K-l 



tdQ(t). 



The average excess holding time for the .4-address case follows from 
(l)as 



Ea= f ld[Q(OY - [ tdQ(i). 



(2) 



The only appntximation in the analysis of the address camp-on model 
is the assumption that the holding time of a line is exponentially dis- 
tributed, so that 



Qit) = 1 -Se 



rft 



(3) 



The holding time is made up of two random variables, namely the mes- 
sage length and the excess holding time. The approximation made in 
(3) implies that the sum of these two random variables is exponentially 
distributed. Tliat this, indeed, is a reasonable assumption is confirmed 






*-z 



VV, A » 



* — z 



W^^A 



ZA,A=dA 



P'ig. 7 — Line seizure sequence. 
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by the close agreement of the delay derived here with the delay derived 
by simulation.' 

Substitution of (3) in (2) gives 

E, = Ab<p r er^' t{l - <pe-'"y'-' dt - <pb \ Ic'" dt 



The above reduces with the binomial expansion f or { I — 3e *" ) to 



For the case of exponential line holding time with mean t and service 
of requests for lines in the order of arrival as is the case here, we must 
substitute in (4) according to Erlang^ 

5^1- Q(0) (6) 

in w^hich c is the number of lines per group, cit the average number of 
requests per line holding time or the total load olTcred per line group, 
and fl/ the average number of requests per message length or the informa- 
tion load. Q(0), the probabihty of no delay, is given by Erlang^ as 

c —an- 

ttr e c 



Q(o) = 1 - — . ,1 °"°:. . ■ m 



are , O-T c 



i_5:«i^ + 



?! c! c — ar 

The unit of time being the average message length renders for (, the 
average line holding time, 

1 + £a = i (8) 

so that 

iJ. ^ ^ - 1. (9) 

a-i 

Substitution of (5), (6) and (9) into (4) brings 

at Oi c — ar ai c — Qt , , 

^^ (A - l\ [Q(0) - If ^' 
'h\ K ) {K + \Y 

With c, ai and A given, we can now compute ar , the total load, itera- 
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tiveiy from (7) and (10) aliove. Then, with Qt known, we find the 
average delay on all messages from (1) and (9) as 

Ui c — Qt ar 

wliich wc recognize as tlie left-hand member of (10). 

It is interesting to take note of the fact that there are generally two 
values of Ot which fulfill (10), Thus, one could conclude that the system 
can operate in two modes, one implying a shorter delay than the other. 
More than one steady state of operation has been observed by others'-^ 
in similar traffic studies. However, it appears questionable that the 
heavy delay mode is stable since the larger of the two Ut which fulfill 
(10) decreases with increasing Ui , which is physically unreasonable. 

APPENDIX B 

Mathematical Analysis of Retrial Model 

When a newly offered message makes its first trial to seize one of c 
lines in a group, let Si ,0 ^ i ^ c, be the vmconditional probability that 
the group is in state i. A lino group is said to be in state * when i out of 
all c lines in the group are busy. The message is reoffered until a line is 
available to each of .1 addressed recei\ing stations; therefore, no mes- 
sages are lost. The load carried on each group equals Ur , the information 
load offered. For the special case of c = 1, we obtain 

-So = 1 - ar (12) 

and 

5] = a,. (13) 

J'or c > I, i.e., for more than one line per group, .S', depends not only 
on a, , hut also on the procedure by which lines are made busy. By 
procedvu-e is meant the type of distribution of the length of the intervals 
between fine re(iuests, and whether unsuccessful attempts form queues 
or are withdrawn. For Poisson input at the rate of Oo and withdrawal of 
blocked attempts, Erlang loss probability^ gives 

S, = -^^^^ (14) 



in which 



a. = j— ^ do) 
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In the switching model considered here, the total input to each line 
group, i.e., the total load offered, is not Poisson distributed and its mag- 
nitude differs from ao . This is so because repeated attempts are blocked 
with a probabiUty which is larger than Sc and because for A > 1 blocking 
to any one of the A addressed receiving stations will cause a retrial. In 
order to compute Si , it is assumed that for a sufficiently large retrial 
interval r, the total input will at least resemble ao ui distribution and 
magnitude. This approximation may be justified for Hght Hne occu- 
pancies but it becomes increasingly unrealistic with increasing line oc- 
cupancies. The approximation for ao , therefore, is used here only to com- 
pute values for Si when c ^ 2. For c = 1 the values for Si are exactly 
determined by (12) and (13) above. For c > 1 we solve (14) and (15) 
iteratively with i = do obtain ao and then solve (14) to obtain approxi- 
mations for Si when c ^ 2. The approximations for S, so obtained are 
used for the unconditional state probability both in the one and three- 
address cases. 

Next, we will consider conditional probabilities which take into ac- 
count past known states of a line group. Let it be known that at a given 
time (o there are j lines busy in a group, ^ j ^ c; what then is the 
probability that at U -{- t there are i lines busy? This conditional prob- 
abihty is called Xij(t). For Poisson input and exponential line holding 
time, Xi,j{t) is given by a well-known set of first-order differential equa- 
tions.* Here, now, we must take into account that the superposition of 
first and repeated attempts results in an input which is not Poisson. 
We let w(t), /o ^ ^ ^ 'o + T, be the instantaneous line request rate or 
the density of requests. As was said above, X,j{t) is conditioned on 
state j of the group at h . Consequently w(0 depends also on the state 
of the group at /o , and this important point should be kept in mind, 
particularly since the notation does not always remind the reader of 
this condition. 

Assume for the time being w(l) is known for every value of ( in the 
interval (/o , 'o + r). The differential equations defining X,,,(/), to ^ 

i ^ ^ + T, are 

Xo,/(i) = -o,(t)-XoAt) '\- x,M) 

Xi./it) =co(0-X._i,KO- [«+a;(0]-X-.;(O + 



(16) 

(i-F-l)-X,+i,,(0 



f or < I < c 

X,,/{t) = oi(t)-X,-ij(t) - c-X.i(0- 
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The condition that j lines are busy at time Iq is taken into account by 

1 1 for i = j 

XiM«) = { (17) 

[0 for i 9^ j, 

and, as said before, by w(/). For the special case i = U -\- t, the argu- 
ment is dropped and X.jiU + t) is abbreviated to Xij . The system 
(16) can be solved with Laplace transforms for c ^ 2, and by numerical 
integration for c ^ 8. 

Two cases are considered in the following: one in which each message 
has one address, i.e., the case .4^1; the other in which each message 
has three addresses, i.e., the case A = 3. 

B.l One Address per Message 

In order to compute delay for the case A = 1, we must know the 
probability of finding the line group in state c at /o + t, given a state c 
at (q jIq — T, to — 2t, etc. In other words, we must know the probabiUty 
that a message is blocked twice, three times, four times, etc. To simphfy 
the notations for the case /I = 1, we write C, for the conditional prob- 
abihty that all c lines of a group are busy at U + r, given a state c at fo , 
ta — T, • ■ ■ , to — ir. Ci , therefore, denotes the probability that a one- 
address message is blocked i 4- 2 times in a row. Co is identical with 
Xc.c and can be computed by (16) provided w(0 is known. It will help 
to keep matters clear if, for the case A ^ \, w(0 is subscripted so that 
o}i{t) refers to the condition that at to , to — r, ■ ■ ■ , /o — ir all c lines of 
the group are known to be busy. For instance, the line request rate used 
in (16) to compute Co is called «o(0- Values for C, , t ^ 1, are computed 
from (16) in the same manner as X^.c , except that w,(i) is conditioned 
as indicated later. Hence, the numerical values for Xc.c obtained from 
(16) with uii(l), i ^ 0, are equal to C,- . 

Let us now discuss the procedure by which wi(/) is obtained for the 
case A = I. We will find functions Li which are conditioned on a state 
c Sit to , to — T, ■ ■ • , /o — ir, such that 

a.,(0 = L,(ro,ri,c,, ---.o. (18) 

This means that we cannot obtain an exphcit expression for w,(0 
since aj,{0 is needed to compute C, from (16). But with (18) we come 
into a position which allows us to assume values for C, , compute u}i(t) 
from (18) and then use the so-computed w,{() to obtain Ci from (16). 
Through successive iterations stable solutions are obtained for C, such 
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that (19) is satisfied 

^„ < Co < C'l < G < ■■• < 1. (19) 

The above approach, of course, is extremely tedious when C. must be 
computed for large numbers of i. The problem is simplified when one 
assumes for large enough i that d = r,+i . Computations have shown 
that for retrial intervals t ^ 0.5 we may reasonably approximate d = 
C,+i when t S 4. This approximation, consequently, is used in the 
one-address case. It is mainly for this approximation that the analysis 
for the case A = 1 is Hmited to retrial intervals t ^ 0.5. 

Before we define the function Li in (18), we will give the method by 
wliich the successive iterations are performed by a computer program to 
compute C, . We start by iterating for a stable value of C\ with d = Co 
for i ^ 1. Next, we iterate for a stable value of Ci with Co fixed and 
d = Ci for i ^ 2. Now we go back and iterate for a new value of Co 
with Ci fixed and C, - Ci for i g 2. This process is continued until no 
further changes in Co and Ci are detected. Continuing one step further, 
we iterate for a stable value of Co with Co and Ci fixed and d = C. for 
1^3. Again, we back up and search for a new value of Co with Ci and 
Ca fixed, then search for a new value of d with Co and Ca fixed and finally 
search for a new value of d with Co and Ci fixed, all with d = Ca for 
i ^ 3. We proceed in steps in tl^e manner described above until finally 
no changes are detected in Co , Cj , C. , C, , d with d = C^ for i § 5. 

With the approximation d = C^ for i ^ 5, we can write for (18) 

c.,(0 = L,(Co,C:,C2,C3,C4,0- (20) 

Line requests are made by first and repeated attempts. First attempts 
arrive independent of time with a density );. Repeated attempts arrive 
with a density iii(t), in which i refers to the condition that all lines arc 
busy at time l» , /<, - t, ■ ■ ■ , /o - «> and / is some time such tliat 

to ^ i ^ lo -\- r. With these definitions we substitute for L,-(Co , Ci , C2 , 
C3 , C4 , /) in (20) 

aj.(0 = V + Ui{t). (21) 

The density of first attempts, according to definition, is numerically 
equal to the information load offered or 



V = Ui . 



The density of repeated attempts in the interval (^o , ^o + t) is derived 
from first attempts which are made before to and are blocked. For in- 



(22) 
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stance a kth attempt, k ^ 2, occurs at t, to ^ t ^ to -{- t, if the attempt 
occurs first at/ — (k — \)t and all lines are busy at ( — {k — 1 ) r, 
t~{k—2)T, ■ ■ ■ , (— T. As an abbreviation we writePr (I, |/i,(2 , ^3 , ■ ■ ■ ) 
for the probability that all lines are busy at t^ , conditioned on all lines 
busy at (1 and (2 and tz • ■ • -As before, let t be an instant in time such 
that to :^ i ^ to -\- T. With the condition that all fines are busy at to , 
to — r, ■ ■ • , to — IT, we obtain for the density of repeated attempts at t 

Ui(t) = aiiVr {t — T \ to , ta — T, ■ ■ ■ , to — ir) 

+ Vv {i - 2T\to,to- T, ■■■ ,to- ir) 

•Vt {i - T\to,to - T, ■■■ ,to- ir, t - 2t) 

+ Pr (( - 3r |/o,/o - T, ■■■ , /o - ir) 

■Ft (t - 2T\to,fo- T, ■■■ ,to- IT, t - 3t) 

■Pr (i - T I /o , ^n - T, ■ • ■ , (0 - IT, ( - 3t, i - 2t) + • ■ •] . 

We are left with the problem of expressing Vr {t^\ ti , t2 , ts , ■ ■ •) m 
the above as functions of Co , Ci ,€2 , C3 and C* . Assume that symmetry 
exists such that for any positive length of time I 

Pr (t,^l\t,,t,-T,---,t,- kr) 

= Ft {t,- l\t, ,(;-{- T, ■■■ ,t,-\- kr). 

In the above it is assumed that traffic congestion builds up to an afi- 
fines-busy condition at ti . I^ -{- t, ■■-, t^ -\- kr in the same manner as it 
subsides after tj; + kr. This assimiption may not be exact for the retrial 
system but this concept is used here since it is expected to give a good 
enough approximation for the folIo\\ing reason. 

If a group is bvisy, say, at ta , then it must be expected that part of the 
traffic which (contributes to the congestion at to is reofTered traffic. The 
fact that congestion occurs at to implies that all lines were busy at 
to — T, to — 2r, etc., with a larger probability than indicated by the un- 
conditional state probal>ility Sc . As an approximation to the function 
by which traffic is expected to build up we construct linear functions 
in time. For example, we assume that the probability of blocking at 
some time l^ < l,i builds up to an all-lines-ljusy condition at l» — t and 
to as shown in l^g. 8. Also, we assume independence of events that are 
not really independent. lor instance, we assume that blocking between 
to — T and /(I occurs with a probability lFi.i(0 as defined below. Similarly, 
independence is assvuned between events occurring witli probability 
N,(l), Mi(t) 01- Wi.jd) and the event which causes a repeated attempt 
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Fig. 8 — Sample of build-up function. 

at some time prior to to . The approximations made in the expressions 
below may account for some of the differences which are observed be- 
tween theory and simulation. 

A first attempt, made at some time / — jr, In ^ I ^ io + r,j ^ 1, is 
blocked in the nomenclature of (22) witli a probability 

Pr (/ - JT \ to , t-a - T, ■ ■ ■ , lo - ir) 

for which we approximate 

Ni(l) = - U + C\ir - i)] for j = i + 1,1^0 

T 

Miit) = - [CV + SAt - t)\ for; = i + 2, i ^ 

T 

and with Ni*(t) = (\/rMCi - 1) + t] 

Wu(t) = 1 - [1 - Nj-r(mi - Ni-j*{t)\ 

for 1 ^ j ^ i, i ^ 1. 

A fcth repeated attempt made at some time t — jr, to ^ i ^ to -\- t, 
j ^ 1, is blocked in the nomenclature of (22) with a probabihty 

Pr [i - JT\to,tt, - T, ■■■ , to - IT, 

t- (j+ 1)t, i - (j + 2)t, ■■■ ,t - U-^ ^)t1 
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for which we approximate 

1 - (1 - C._i) [1 - Ni(t)] fori = i+l,i^O 

1 - (1 - C,_i) [I - M,(0] fori = i + 2,i^0 

Ck-i forj ^ I + 3, z ^ 
and finally 

1 - (1 - Ck-i) [1 - Wij{t)] for 1 ^j^i,i^l. 

We make use of the above expressions as shown in (22) to obtain 
Ui(i). The subscript i of U;(t) corresponds with tlio value of i in the 
above approximations and refers to the condition that all Unes are busy 
bX to , to — T, ■ ■ ■ , lo — ir. For every such i, i = 1, 2, 3 and 4, we have, 
according to (21), a L,(Co , Ct , C2 , Ca , C4, t) and an t.J,(/), and can per- 
form the iterations outlined before to compute Co , f'l , €■> , C3 , and C4 . 

Continuing in the analysis of the one-address case, we will now 
evaluate the delay. The probability that a one-address message is de- 
layed exactly ir is given by Di{i) which is 

AfO) - 1 - .S, 

Difl) = S,{\ - Co' 

i— 2 

D,{i) = Sr(l - C,_i) riCr i ^ 2. 

Hi{i), the probabihty that the delay is greater than ir for the one- 
address case, is 

Hi(i) = 1 - T.D,(k), i^ 

k=0 

which reduces to 

i/i(0) = S, 

Hi{i) = Sc ncj, i^ 1. 

The average delay is obtained as the summation of all possible delay 
values multiplied by their respective probabiUty of occurrence and is 
given by 

di = T ^i Di(i) 
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which reduces to 

Ch- rS,\l+ Z nCy]. (23) 

The above formula is used to compute the average delay on all messages 
for the case A = I. The unconditional probabihty of finding all lines 
busy, represented by Se , and the conditional blocking probabilities d 
are approximated liy the methods outlined before. 

R.2 Three Addresses per Message 

The delay for the case A = 3 is computed in a Markov process. This 
means that we are considering only a lirst-order dependency, since we 
make the assumption that the conditional probability Xij of finding i 
lines of a group busy at io + t depends only on state j of that group at 
/„ . This and the following approximations appear justifiable in the multi- 
address case when considering the multitude of factors which determine 
the hne request rate between (o and to + t. The principal assumption 
for the case ^ = 3 is that oj, the sum of the densities of line requests 
of first and repeated attempts during /o and /„ + r, is nearly Poisson dis- 
tributed and therefore independent of any states at or before h> . llecall 
that fov A = 1 we have been concerned only with the conditional prob- 
alnlity of state c at to -\- t given also u state c at to or at to and to — r, 
etc. For the case A = 3, however, we are concerned with an .Y,,> for all 
values of i and j, O^i^c, OSJ^ c, as will become apparent later. 
For a known o:, we obtain Xij from (10). The condition that j trunks 
are busy, now, is accounted for only by the initial condition as given in 
(17). The density of line requests u is obtained similarly to (21) as the 
sum of the densities of first attempts ar , and of repeated attempts u 
which, according to our assumptions for A = d, are time independent. 
It is obvious that in the case of A addresses per message, A > \, a hne 
request in a given group is made only when the condition is fulfilled 
that the remaining A - 1 groups are not busy. Since mdependence is 
assumed wc can .set (1 — S^Y for this condition in the three-address 
case and obtain 

a: = (a, + w)(1 - S.f (24) 

for the density of line requests in the interval (^o , fo + t) in any given 
line group. The expression given in (24) above, of course, is an approxi- 
mation since in reality u is dependent on the state of the group at (o 
and since independence is assumed between the event causing a repeated 
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attempt during (to , to -\- t) and the event causing all lines to be busy 
in the other two addressed line groups. But, as said before, these de- 
pendencies are believed to be noncritical for the multiaddress case, so 
that w is considered to be independent of time. 

As in the one-address case, we are left with the problem of defining w, 
the density of repeated attempts, which is expressed below as a fvuiction 
of Xe.c . We obtain u by the following approach. Let G bo the probabihty 
that a first attempt Is blocked at some time prior to /c . This probabihty 
is approximated by 

a - 3.S,(1 - SS' + 3.S/(1 - ,S,) + ,S,'. (25) 

The above is an approximation be<iause it assumes independence be- 
tween the event which causes a group to be in state./, ^ ;/' ^ c, at to 
and the event causing all lines of a group to be busy at some time prior 
to /i) . For the probabihty that a Wn. attempt, k g 2, is blocked prior to 
to we approximate 

R = A%.r(l - SS- + i^.S.d - A',..)(l - N,) 

+ 2A',.,-SV(1 - .S',/) + .S',-{ I - A,,,) + A,,,,SV 

for which it is assumed that at the /; — 1st attempt one hue group was 
in state c, i.e., busy, but without having made any assumptions about 
the state of the remaining two groups. The expression given for R in 
(20) is an approximation .since, as before, the known state of a group at 
tft is ignored and since only a first-order dependency is considered, as 
mentioned earlier. The flensity of repeated attempts is obtained simi- 
larly to the one-address case by considering all attempts which were 
blocked prior to to so that 

u = a, a -\- (i,GR + ajGR- + ■ ■ • 

or 

Substitution of (27) in (24) gives 

R in (28) abo\-e is a function of A^.^ . This means that we cannot find w 
explicitly since co is needed to compute A',., as outhned in (16). To find 
w we again must iterate by assuming a value for X^.c in (26), recompute 
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Xc.c from (16) and then use the recomputed X^.c in (26). After having 
found a stable value for o: we can, from (16), readily compute Xij for all 
^ I ^ c, ^ J g c. 

In order to compute the delay for the three-address case we con- 
sider Pi,j.k{k), which is defined as the probability of finding the 
three groups in state i, j and h respectively at the k + 1st attempt, 
k ^ 0. Pij,h(k) is obtained recursively by finding ah possible ways in 
which the states of the three groups have changed to states i, j and h at 
the k + Lst attempt when at least one group was busy at the /cth attempt. 
Using the approximation of a first-order dependency, as mentioned 
before, we get for fc ^ 1 

PijAk) = i: J2Pc,rAk - 1)-A%vZ,>-Xm 

+ E EP..c,„(/o - l)-X,,,-XiyXn, (29) 

+ ILEP'i.rAk - l)-Xi,,-Xi.r-X,.,. 

3=0 r=0 

At the first attempt, 

Pi.iAQ) = SSjS, . 

A tlu-ce-address message succeeds at the kih attempt, k ^ 0, when at 
the fcth attempt all three groups are in states other than c. The prob- 
ability of a delay of exactly kr, then, is given for the three-address case 

by 

1=0 3=0 ft=-0 

Hz{k), the probability that the delay is greater than kr for the three- 
address case, is 

Hz{k) = 1 - HB.ii), k^O. 

i=0 

As in the one-address case, we find the average delay on all messages for 
the three-address case by summing over all possible delay values multi- 
plied by their respective probability of occurrence 

00 

^3 = rT,k-D,(k) 
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or with (30) 

(h = rfl Z Z Z lc-Pij.>.(k). (31) 

A=i 1=0 y=o /i=o 

Recall that Pij.hik) is obtamed recursively as shown in (29). The un- 
conditional state probability S, and the conditional probability Xij 
which are both needed in (29) are approximated as described earlier. 

APPENDIX c 

Mathematical Analysis of Message Storage Model 

The delay for the message storage model is computed by well-known 
methods of traffic theory and is given here only for reasons of com- 
pleteness. The delay in the delivery of a copy of the message to a given 
station is — according to the switching model — independent of the 
delay in the delivery to any other station. Delayed messages form queues 
in the order of arrival. An analysis for queued service and exponential 
Hne holding time was made by A. K. Eriang. 

According to Eriang^ we find for the average delay on all messages to 
any given receiving station 



d - F(0) (32) 

c — ai 



with 



no) 




c —ar 

tti e c 
c\ c — a, 




1 - 


^ fl/'f""' ajV' 


c 



i! c! c — ai 

The delay distribution, expressed as the probability that the delay is 
greater than /, is computed from 

The curves for the message storage model are calculated from (32). One 
should, however, bear in mind that in certain specialized applications of 
data communication a copy of the message must sometimes have been 
delivered to all addressed receiving stations before the message is of use 
to any one station. One would then be interested in the average delay 
until a line is found to the receiving station with the longest delay of all 
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Stations addressed by the message. For A addresses per message this 
delay is given by 



'0 

which for -1 = 3 reduces to 



d^ ^ [ id[l - F(t)V, 



„.Z(0li3-3^(0)+M.. 
c — a, [ 2 J J 
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